Many-body physics in the classical-field description of a degenerate Bose gas 
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The classical-field formalism has been widely applied in the calculation of normal correlation functions, 
and the characterization of condensation, in finite-temperature Bose gases. Here we discuss the extension of 
this method to the calculation of more general correlations, including the so-called anomalous correlations 
of the field, without recourse to symmetry-breaking assumptions. Our method is based on the introduction 
of U(l)-symmetric classical-field variables analogous to the modified quantum ladder operators of number- 
conserving approaches to the degenerate Bose gas, and allows us to rigorously quantify the anomalous and non- 
Gaussian character of the field fluctuations. We compare our results for anomalous correlation functions with the 
predictions of mean-field theories, and demonstrate that the nonlinear classical-field dynamics incorporate a full 
description of many-body processes which modify the effective mean-field potentials experienced by condensate 
and noncondensate atoms. We discuss the role of these processes in shaping the condensate mode, and thereby 
demonstrate the consistency of the Penrose-Onsager definition of the condensate orbital in the classical-field 
equilibrium. We consider the contribution of various noncondensate-field correlations to the overall suppression 
of density fluctuations and interactions in the field, and demonstrate the distinct roles of phase and density 
fluctuations in the transition of the field to the normal phase. 

PACS numbers: 03.75.Hh 



I. INTRODUCTION 

The experimental realization of Bose-Einstein condensa- 
tion (BEC) in dilute atomic gases (THS has lead to a resur- 
gence of interest in the theory of weakly interacting Bose 
gases [4-6 1. Theories of weakly-interacting BEC were ini- 
tially developed [7, 8 1 in the hope of obtaining insight into the 
physics of the strongly interacting superfluid helium, and can 
only give a qualitative account of the physics of the liquid su- 
perfluid phase (5). By contrast, the advent of dilute, weakly 
interacting Bose gases in the laboratory provides for the di- 
rect comparison of theories of weakly interacting BEC with 
experiments (for a review see Ref. (6)). More generally, these 
systems offer an unprecedented opportunity for the quantita- 
tive experimental evaluation of quantum-field models at finite 
temperatures, and away from equilibrium, as noted by other 
authors flOHBI . 

The analysis of the weakly interacting Bose gas at finite 
temperature is significantly complicated by the necessity of 
considering interactions between excitations of the conden- 
sate. A self-consistent mean-field approach, based on ap- 
proximate factorizations of field-operator products [14| leads 
to the so-called Hartree-Fock-Bogoliubov (HFB) formalism 
for the finite-temperature field, which is also obtained in a 
variational approach to the problem |[T5l [161 . This descrip- 
tion, however, violates known exact constraints on the exci- 
tation spectrum |17|, due to the inconsistent introduction of 
(leading-order) many-body effects lfl8l into the description of 
atomic interactions. The simplest resolution of this problem 
is provided by the so-called Popov approximation lfT4l [T9l to 
the HFB theory (HFB-Popov), in which interactions between 
excitations at the Hartree-Fock level are retained |4|, while 
many-body corrections are neglected in the treatment of all 
interatomic interactions. 



A systematic perturbative treatment of the beyond- 
quadratic terms of the Bose-field Hamiltonian [20 1 (see also 
Refs. |01|2Tj) demonstrates that the failure of the HFB method 
results from an inconsistent treatment of Hamiltonian terms 
cubic in the Bose field operator in the mean-field factoriza- 
tion approximation. This provides support for an extension of 
the conventional HFB approach [22 23 1 in which the Popov 
theory is augmented by spatially dependent effective interac- 
tion strengths, which serve to upgrade the Hartree-Fock inter- 
actions between atoms to interactions mediated by a many- 
body T matrix, which accounts for all ladder-diagram pro- 
cesses contributing to interactions in the finite-temperature 
medium (4] [TS] |24l l25ll . However, in this theory the zero- 
energy limit of the T matrix is substituted for the interac- 
tion strength in all condensate-excitation interactions, neglect- 
ing the dependence of the T matrix on the collision ener- 
gies. More complicated (bubble-diagram [26|) effects which 
appear at the same order of perturbation theory |20| as the 
ladder-diagram corrections are also neglected. Kinetic the- 
ories which include processes beyond those contained in the 
HFB theories, such as the exchange of atoms between the con- 
densate and thermal cloud, have been considered by several 
authors lfl6l I27H321 . 

An alternative approach to modeling the finite-temperature 
Bose gas is provided by the so-called classical-field (or c- 
field) formalism [33-36]. In this approach, one considers 
the dynamics of a classical dynamical system correspond- 
ing to the high mode-occupation limit of the formal second- 
quantized field theory lfT5l |37| for the dilute Bose gas. The 
classical-field model provides a leading-order description of 
the long-wavelength physics which dominate the critical be- 
havior associated with the Bose-condensation phase transi- 
tion 11381 . More generally, the dynamical classical-field equa- 
tions of motion arise as the "classical" component of the 
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atomic-field evolution in the Wigner representation of the 
second-quantized field [39-42]. The classical-field model is 
thus expected to give a good description of the most highly 
occupied modes of the system, including the condensate and 
its low-lying excitations, for which quantum fluctuations can 
be safely neglected. The introduction of an explicit single- 
particle energy cutoff restricts the system to these low-energy 
modes Il43l [44). In closely related stochastic-field meth- 
ods Hp 45 ], the description of the low-energy modes includes 
explicit damping and noise terms arising from the coupling of 
these modes to the eliminated high-energy component of the 
Bose field. The classical wave-mixing dynamics of the result- 
ing (Hamiltonian or stochastic) field equation of motion then 
provide an intrinsic nonperturbative description of all inter- 
action processes in the low-energy field, neglecting only the 
quantal nature of excitations, which is essentially irrelevant in 
the high-temperature regime. The inclusion of many-body ef- 
fects beyond Bogoliubov theory in the Hamiltonian classical- 
field theory, and agreement of this theory with the second- 
order perturbative treatment of Ref. GUI was demonstrated 
for a homogeneous field in Refs. Il43ll46l . 

In this article, we present a comprehensive, explicit demon- 
stration that the equilibrium classical-field dynamics pro- 
vide an intrinsic description of many-body interactions in the 
finite-temperature, harmonically confined Bose field. Intro- 
ducing appropriate U(l)-symmetric field variables, analogous 
to the modified ladder operators of number-conserving Bo- 
goliubov theories [47— 50|, we quantify the anomalous and 
non-Gaussian nature of the field fluctuations. These corre- 
lations reveal signatures of many-body processes neglected in 
(e.g.) the HFB theories. We discuss the importance of these 
processes in shaping the condensate orbital, and demonstrate 
the consistency of the Penrose-Onsager definition of conden- 
sation in the classical-field equilibrium. In particular, we show 
that the Penrose-Onsager condensate appears as an effective 
eigenfunction of the total mean-field potential it experiences, 
provided that the contributions of anomalous averages (both 
pair and triplet) to the latter are taken into account. We also 
consider the temperature dependence of the many-body ef- 
fects, and discuss the relation between condensation and the 
overall suppression of density fluctuations in the system. Our 
results reveal that the partially condensed Bose gas exhibits 
non-trivial correlations as a result of purely thermal (classi- 
cal) fluctuations of the atomic field, which dominate over the 
quantum fluctuations of the field in realistic experimental sys- 
tems eh. 

This article is organized as follows: In Sec.|ll]we briefly re- 
view the derivation of the projected Gross-Pitaevskii equation 
(PGPE) formalism, and give the parameters of the system we 
investigate. In Sec. Ill we describe the microcanonical inter- 



pretation of the PGPE, and define U(l)-symmetric classical- 
field variables appropriate for the evaluation of anomalous 
correlation functions of the classical field. In Sec.[lV]we dis- 
cuss the local correlation functions of the noncondensed com- 
ponent of the field and their relation to many-body interac- 
tion effects described by the classical-field model. In Sec. [V] 
we consider the dependence of the noncondensate correlations 
and interaction effects on the energy (and thus temperature) of 



the classical-field equilibrium. In Sec. VI we consider the var- 
ious contributions to the suppression of density fluctuations 
in the field which are neglected in mean-field theories. In 



Sec. VII we summarize and present our conclusions. 



II. FORMALISM 



A. PGPE formalism 



A detailed review of the formalism of (projected) classical- 
field methods has recently been given in Ref. [33 ], but for the 
reader's convenience, we briefly describe the relevant details 
of the formalism here. 

Formally, the physics of the harmonically trapped dilute 
Bose gas is governed by the second-quantized Hamiltonian 

+ ^JdxJ rfx'4' t (x)4' 1 (x')t/(x-x')4'(x')4'(x), (1) 

where the single-particle Hamiltonian is 

-H 2 V 2 m 



H, 



sp 
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[u 2 x x 



2 + u 2 y 2 + oj 2 z 2 ], (2) 



and t/(x) is the exact interatomic potential. We introduce a 
single-particle subspace L spanned by eigenmodes T„(x) of 
the single-particle Hamiltonian [H sp Y„(x) = e„Y„(x)] with 
energies e„ less than some cutoff E m . dx , and a complemen- 
tary subspace comprised of the remaining high-energy modes. 
Provided E max is chosen such that the high-energy modes are 
essentially unoccupied, the dynamics of these modes can be 
integrated out to obtain an effective Hamiltonian for the low- 
energy (coarse-grained) Bose field *Tl(x) = ZneL <2 n T„(x), as 
shown by Morgan iTZOl . Atomic interactions described by the 
effective Hamiltonian are mediated by an approximate 2-body 
T matrix, and the interaction can thus be rigorously approxi- 
mated by a "contact" potential, with a renormalized coupling 
constant C/q. The low-energy Hamiltonian then takes the form 



fdx%(x)H sp %(x) 

+ Vl J" dx¥ L (x)¥ L (x)y L (x)V L (x), 



(3) 



which defines an effective field theory [52| for the coarse- 
grained field T'iXx). 

We then further divide the low-energy region L into a co- 
herent region (or condensate band) C = {« : e„ < e cut }, 
spanned by single-particle eigenmodes Y n (x) with energies be- 
low some classical-field cutoff e cut , and a complementary in- 
coherent region I = {n : e cut < e„ < E max }. Introducing the 
projector 

p{f{x)} = ^ Y n (x) f dy T„*(y)/(y), (4) 
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onto the coherent region C, we define a coherent-region field 
operator 



ifcx)=P{%(x)} = Y j a n Y n (x). 



(5) 



neC 



Neglecting the coupling of this field operator to modes of the 
field in the incoherent region I, we find that ^i{x) is gov- 
erned by a Hamiltonian obtained from Eq. Q by the re- 
placement T'l(x) — > iA(x). The classical-field approxima- 
tion is then made by demoting the operators a n in Eq. (|5]l 
to classical variables a n , thereby defining the classical field 
|A( X ) — Znec <XnYn(x)- The evolution of the field \f/{x) is gov- 
erned by the classical-field Hamiltonian 



HcfW = Jdxf(x)H sp i(,(x) + 



u -i\m\\ 



(6) 



and Hamilton's equations for the classical variables a n ob- 
tained from Eq. (|6]i can be expressed concisely as the field 
equation of motion 



in ^T = 9 {( Hs p + ^ol^Cx)! 2 ) «A(x)} , 



dt 



(7) 



which is termed the projected Gross-Pitaevskii equation Il33l . 
Equation |7]i describes the classical (thermal) fluctuations of 
the coherent region, while neglecting the contribution of quan- 
tum (vacuum) fluctuations, which is a valid approximation in 
regimes of significant thermal excitation, where mode occupa- 
tions are large, and thermal fluctuations dominate. For an ap- 
propriately chosen cutoff e cut , one can augment the classical- 
field description of the coherent region C with a mean-field 
(Hartree-Fock) description of the incoherent region I, allow- 
ing for accurate quantitative predictions for the full Bose-field 
system [53]. In this article, we consider only the correlations 
of the coherent region C that arise from the kinetics of Eq. (FT). 



B. System parameters 

For the purpose of the numerical implementation of the 
PGPE, we rescale Eq. (|7]) by introducing units of distance 
ro = yfh/maT r , time to = u>~ 1 and energy eo = hio r , with 
oj r = o) x = a> y the transverse trapping frequency of a cylindri- 
cally symmetric system. With these choices, Eq. (|7]i becomes 



ijL = __vV + -{x 1 +y z + Aiz 1 ^ + 5P{Cn#l», (8) 

where we have chosen the classical field i/r(x) to be unit- 
normalized, and absorbed the total classical-field atom num- 
ber N c into the nonlinear coefficient C n \ - N c Uo/tko r ri. We 

follow Ref. l54l in choosing parameters A, — V8 (a typi- 
cal three-dimensional trap geometry), C n \ — V2 x 500, and 
£cut = 317zo> r . The corresponding ground (Gross-Pitaevskii) 
eigenstate of the system has energy E sa 9N c hu> r . We ana- 
lyze equilibrium configurations of this system with energies 
in the range E e [9.5, 24.Q]N c ha> r . We note that our choice 



of parameters makes no reference to the physical number of 
atoms in an experimental system. The classical-field approx- 
imation becomes asymptotically exact in the "classical" limit 
N c — > oo, U() — > 0, with C„i fixed [33], and our results 
correspond to this idealised classical-limit system. Never- 
theless, for realistic experimental systems in moderate and 
high-temperature regimes, the magnitude of quantum fluctu- 
ations is small compared to that of thermal fluctuations, and 
the PGPE provides a good description of the thermal Bose 
gas 



III. CORRELATION FUNCTIONS 

A. Microcanonical formalism 

A central feature of the PGPE formalism is the micro- 
canonical (ergodic) interpretation of field trajectories at equi- 
librium. It is well known that the solutions of the nonlinear 
Schrodinger (Gross-Pitaevskii) equation exhibit stochastiza- 
tion, leading to the ergodic equilibration of the system and the 
equipartition of the system energy l35l |6TH66l . This prop- 
erty is also shared by trajectories of the PGPE for homoge- 
neous 03]|46]|62) and harmonically trapped H4l 1531 l68l Irj9l 
systems. In the ergodic interpretation, the PGPE trajectories 
provide a sampling of the microcanonical density 



P[*fr,E] = 



const 




HqpW = E 
Ham * E, 



(9) 



defined by the conserved first integrals of the system. As the 
trajectories of i/r(x, t) cover the density P[i[>; E] densely, aver- 
ages in the microcanonical density can be approximated by 
time averages along the trajectories of ifr(x, t). We therefore 
define correlation functions of the classical field as averages 
(or expectation values) of functionals of the field if/(x,t) in 
the density Eq. (|9j, which we approximate by time averages 
of the field (denoted by (• • ■ )). Applying the microcanonical 
thermodynamic formalism of Rugh f70l to the PGPE trajec- 
tories, one finds that the PGPE system evolves over time to 
an equilibrium characterized by a well-defined chemical po- 
tential and temperature [55, 67], providing strong support for 
this ergodic interpretation of the classical-field dynamics. 

A correlation function of particular importance in the appli- 
cation of the PGPE to partially condensed Bose systems is the 
covariance matrix 



G(x,x') = <^(x)^(x')), 



(10) 



which is the classical-field analog of the quantum one-body 
density matrix ll44ll7D . As G(x, x') is Hermitian, we can di- 
agonalize it to obtain a complete basis of eigenvectors l^,(x)} 
(the eigenmodes of the one-body density matrix) with real 
eigenvalues {«,■} (the mean occupations of the eigenmodes at 
equilibrium), such that G(x, x') = Tii n iX*( x )Xi( x ')- By anal- 
ogy to the Penrose-Onsager definition of condensation IT721 . 
we identify the largest eigenvalue no as the condensate occu- 
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pation, and the corresponding eigenvector xq(x) as the con- 
densate orbital. 



B. U(l)-symmetric correlation functions 

An important feature of classical Hamiltonian systems is 
the relation between symmetries of the Hamiltonian and the 
conservation of quantities during the system's evolution (see, 
e.g., Refs. 1731 1741 ). The Hamiltonian |6]) governing the 
classical-field dynamics is invariant under the U(l) (gauge) 
transformation if/(x) — > 4r(x)e' e , and this symmetry ensures 
that the normalization N[i[f] = Jdx\if/(x)\ 2 of the classical 
field is conserved under the action of Eq. Q. Moreover, as 
Eq. (|6| has no explicit time dependence, the classical-field en- 
ergy HcfW is also conserved under the PGPE evolution. We 
note that microcanonical density [Eq. |9])] of the PGPE sys- 
tem inherits the symmetries of the Hamiltonian (|6j; in partic- 
ular, P[tfr;E] is symmetric under the U(l) gauge transforma- 
tion if/(x) — > if/(x)e' e , and thus only the averages of quantities 
which are invariant under such transformations are nonzero in 
the microcanonical density. 

The U(l) phase symmetry of the classical-field Hamilto- 
nian Hqf is of course also a symmetry of the "fundamen- 
tal" second-quantized Hamiltonian ([TJ, in which case it cor- 
responds to the conservation of particle number under the 
action of the corresponding Heisenberg equation of motion 
for 'i'(x) 1151 . However, in traditional mean-field theories 
of Bose condensation, the condensate is assumed to acquire 
a definite phase, breaking this symmetry ll75l l76l . Theo- 
ries built on this assumption do not strictly conserve parti- 
cle number, and the grand-canonical formalism is typically 
employed to ensure conservation of the mean particle num- 
ber 11771 . Condensation in the field is then associated with the 
appearance of anomalous moments (moments of the field op- 
erator which are strictly zero in a state of fixed particle num- 
ber) such as the condensate mean-field (¥), and the anoma- 
lous thermal density (66) (where 6 = 4* - ("¥)). In Ref. l54l . 
it was demonstrated that although the U(l) symmetry of the 
Hamiltonian Hqf is formally inherited by the microcanoni- 
cal density Eq. (|9j, the classical-field solutions do exhibit a 
symmetry-breaking aspect, which allows for the calculation 
of anomalous averages such as (iff) and (66) in analogy to 
symmetry-breaking approximations to the second-quantized 
field theory. 

An alternative approach to the theory of Bose condensates 
that respects the U(l) symmetry of the quantum-field Hamil- 
tonian, and therefore conserves the number of atoms in the 
system, was presented by Girardeau and Arnowitt B71 . and 
later rediscovered by Gardiner [ 48 1 and Castin and Dum 1 49 1 . 
This approach is distinguished from the symmetry-breaking 
approaches in that no spurious phase is assumed for the con- 
densate, and the fluctuations of the quantum field about the 
condensed mode are described in terms of modified ladder op- 
erators hi ~ [ag/(aQfl )^ 2 ]a/ (and their Hermitian conjugates) 
f78l . which demote a field quantum from a noncondensate 
mode <pj(x) into the condensate mode (f>o(x) (and vice versa). 
In effect, these operators serve to describe the excitations of 



the system in a picture in which the (indeterminate) phase of 
the condensate mode is canceled, providing a rigorous basis 
for the construction of a theory of the fluctuations around the 
condensate. Such an approach has been used by several au- 
thors in developing theories of Bose-Einstein condensation at 
zero and finite temperature 



In this article, we take an analogous approach to character- 
izing fluctuations of the classical field about the condensate 
mode. We introduce the fluctuation (or noncondensate) field 



A(x, t) = 



«n(0 



y/a* (t)a (t) 



6iff(x, t), 



(11) 



where ao(t) = Jdxx* (x)iff(x, t) is the classical-field amplitude 
corresponding to the condensate mode ^o(x), and 

6iff(x, t) = iff(x, t) - Xo (x) J dx'xoix'M*', t), (12) 

is the component of the classical field orthogonal to the con- 
densate mode. In terms of the Poisson brackets (73] |74) de- 
fined as 



{F,G} 



,j d x 



6F 6G 



6F 6G 



6iff(x) 6iff* (x) 6iff* (x) 6iff(x) 



(13) 



where the projected functional derivative operator 6/6if/(x) 
Unec Y*(x)d/da n 1331 . we easily find 



{A(x),AW} = 0. 



(14) 



The field A(x) is therefore formally invariant under global 
rotations of the classical-field phase. Thus, whereas mo- 
ments of the noncondensate field which are not invariant un- 
der such rotations (e.g., (6ift6iff)) necessarily vanish in the mi- 
crocanonical ensemble due to the U(l) phase symmetry, the 
analogous moments of A(x) (e.g., (AA)), which we will re- 
fer to as anomalous moments of A(x) in the remainder of 
this article, may legitimately acquire nonzero values in the 
ensemble. In the limit of a perfectly coherent condensate, 
for which (Iffol 2 ") = (Icol 2 )"* it is clear that the first moment 
(A) = (a* 6if/)/ yfhj) = 0, as is appropriate for a fluctuation 
variable. It should be noted that more generally the vanishing 
of (A(x)) is not guaranteed a priori [79|; nevertheless we find 
(A) = to good accuracy in our simulations. We note that 
normal correlation functions [those which are not anomalous; 
e.g., <|A(x)| 2 )] are manifestly real, i.e., it is clear a priori that 
complex conjugation has no effect on such correlation func- 
tions. Anomalous correlation functions, by contrast, are not 
fundamentally constrained to be real in this manner. However, 
we find that the anomalous correlation functions we consider 
in this article are predominantly real quantities (relative to a 
real condensate orbital), as a result of detailed balance in the 
equilibrium field. 
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1. Numerical procedure 

Our procedure for calculating correlation functions is as 
follows: We form random initial states B3l l44l with pre- 
scribed classical-field energies. After evolving these initial 
states to equilibrium, we construct the one-body density ma- 



from the critical regime associated with the transition to the 
normal state 



trix [Eq. ( 10 1] by ergodic averaging of the field trajectory. The 



bulk features and simplest correlations of the field equilibrate 
on a time scale of ~ lOOwT 1 ||331|54 , 68 1; however, we do not 
begin our ergodic averaging of the field trajectory until after 
an equilibration period of 1 200<w7 1 , to ensure that all moments 
of the field have settled down to their equilibrium values. We 
then form ergodic averages from 1 .8 x 1 4 equally spaced sam- 
ples of the classical field taken over a subsequent period of 
7200wr' of the field evolution. We diagonalize G(x,x') to find 
the condensate orbital XoOO an d its mean occupation no. For 
the irrotational system we consider here, the condensate or- 
bital has (aside from small numerical fluctuations) a uniform 
phase. Before proceeding we absorb the overall phase of^o(x) 
into a global (time-independent) phase shift of the field tra- 
jectory used to construct the microcanonical ensemble. This 
amounts to making a choice of gauge such that the condensate 
orbital is real and positive, and in terms of this orbital we de- 
fine the condensate wave function <I>o(x) = y/rioxoi^)- From 
each member of the ensemble we then form A(x) according to 



Eq. (Ill, and we average products of A(x) over the ensemble 



to form its moments. We note that Eq. ( 1 1 1 is singular when 
the overlap of the classical field ifr(x, t) with the condensate or- 
bital ^o(x) vanishes. However, the configurations of the field 
for which ao — are a set of measure zero in the microcanoni- 
cal density, and in practical simulations we always find a finite 
value for ao, even in the high-temperature regime where the 
mode^o(x) is incoherent and undergoes large number fluctu- 
ations [80|. 



IV. FLUCTUATION CORRELATIONS AND 
INTERACTION POTENTIALS 

In this section we characterize the fluctuations of the non- 
condensate field by calculating the local second and third mo- 
ments, and fourth cumulants lIHTI . of the fluctuation field A(x). 
In terms of the U(l)-symmetry preserving approach to cal- 
culating correlations we adopt in this article, we find that 
the field A(x) exhibits fluctuations which are both anoma- 
lous (representing pairing effects induced by the condensate), 
and non-Gaussian (exhibiting deviations from the Gaussian 
ansatz for fluctuations assumed in the HFB theories |[T4l[T5ll ). 
These nontrivial correlations of the field can be related to cor- 
rections to the mean-field potentials experienced by the con- 
densed and noncondensed atoms, due to many-body effects 
fsee e.g.. Refs. 112^1221 1231I2511. 

We focus here on a representative PGPE equilibrium, with 
energy E = 14.5N c hco r , and corresponding condensate-band 
condensate fraction no/N c = 0.50. This equilibrium is there- 
fore a reasonably high-temperature state of the field, well 
above the validity regime of the simple Bogoliubov descrip- 
tion |7J 149] [77J of the noncondensate, while remaining far 



A. Moments and cumulants of the noncondensate field 

1. Second moments 

There are two independent quadratic moments of the non- 
condensate field: the normal covariance matrix (or non- 
condensate density matrix) p(x,x') = <A*(x)A(x')), and 
the anomalous covariance matrix (pair matrix) k(x, x') = 
(A(x)A(x'))- All other quadratic moments can be related to 
these matrices by transposition and complex conjugation. The 
spectral representation of the classical field allows us to calcu- 
late the full off-diagonal structure of these matrices, however, 
in this article we consider only their diagonal elements: the 
normal thermal density 



p(x) = <A*(x)A(x)>, 
and the anomalous thermal density 

k(x) = <A(x)A(x)>. 



(15) 



(16) 



We note that in fact p(x) = (6t//*(x)St//(x)), i.e., this quantity 
can be defined perfectly well without recourse to the consider- 



ations of Sec. IIIB By contrast, the anomalous density a:(x) is 



distinct from the moment (6i//{x)5^(x)), which of course van- 
ishes in the microcanonical ensemble. The physical origin of 
k(x) can be inferred by assuming that a factor 1 /no can be fac- 
tored out of the expectation value Eq. ( 16 1 [which should be 



a very good approximation for a coherent condensate mode 
^fo(x)], whereby we find that k(x) ~ {a* Q a* Q 5il/(x)5il/{x)) / no- 
The anomalous density thus quantifies the correlations of pairs 
of noncondensate atoms with pairs of condensate atoms due 
to the (classical wave-mixing analog of the) Bogoliubov pair- 
promotion process in which two condensate atoms scatter 
each other out of the condensate (and the time-reversed pro- 
cess), which is responsible for the well-known Bogoliubov 
particle-hole structure of excitations in the system ||49l 1771 . 
The anomalous density can also be interpreted as a measure 
of the squeezing of the noncondensate field fluctuations [82]. 

In Fig.[TJa) we plot azimuthally averaged values of the ther- 
mal density p(x) (solid line) and the anomalous thermal den- 
sity k{x) (dashed line) on the z = plane. The condensate 
density |<t>o(x)| 2 (dot-dashed line) is plotted for comparison 
[the peak condensate density |O (0)| 2 = 12.7 x 10~ 3 r Q 3 ]. The 
thermal density reaches its maximum in the periphery of the 
condensate, as is well known from mean-field theories (cf., 
for example, Ref. Il83ll ). The anomalous thermal density k(x) 
is real and negative to within statistical uncertainty due to the 
finite size of the ensemble l54l . as expected, since in equi- 
librium the Bogoliubov pair-promotion process must be bal- 
anced by the corresponding time-reversed process. We note 
that the anomalous density is almost entirely located within 
the extent of the condensate mode, and its magnitude exhibits 
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FIG. 1. (Color online) Correlations of the fluctuation field A(x) 
in a representative PGPE equilibrium with energy E = l4.5N c ha> r : 
(a) Second moments of A(x) (normal and anomalous thermal densi- 
ties), (b) Third moments and (c) fourth cumulants of A(x). 



a small dip around the trap center, in agreement with the re- 
sults of previous studies C2l l23l 1531 l84l . 



the non-Gaussian nature of the fluctuation field A(x). The 
physical origin of the correlator A(x) can be inferred simi- 
larly to that of k(x): assuming a factor 1 / yjh~Q can be fac- 
tored out of the expectation value Eq. (17i, we find that 
A(x) ~ {a* Q 5\jj* Si/jSif/) I yfhj). The function A(x) is thus asso- 
ciated with processes in which two thermal atoms collide and 
one of them is scattered into the condensate (i.e., condensate 
growth processes BT1 1571 [88)). In equilibrium this process 
is balanced by the time-reversed process, which is consis- 
tent with the real value of A(x) we find here, and represents 
the detailed balance of Beliaev and Landau processes (see, 
e.g., Ref. 1 89 1) in the noncondensate, which we emphasize 
are not included in the HFB and HFB -Popov treatments ll20l . 
The correlation function ;f (x) is less straightforward to inter- 
pret ||90ll ; however, analogous correlations are found to enter 
into the equations of motion for the pair matrix and the cor- 
relation function A(x) in generalized mean-field treatments of 
the finite-temperature Bose-gas dynamics lfl6l I3T1 l86l . 



3. Fourth cumulants 

The appearance of the nonzero third moments A(x) and^-(x) 
shows that the fluctuations of the noncondensate field are not 
strictly Gaussian. However, as the classical field we consider 
is weakly interacting, the fluctuations of A(x) should be rea- 
sonably close to Gaussian. We therefore expect the fourth mo- 
ments of the noncondensate field to be approximately equal to 
their naive Gaussian (Wick) factorizations [ 15 1 in terms of the 
second moments of A(x); e.g., 



Third moments 



<|A(x)| 4 > 



2<|A(x)| 2 ) 2 + |<A 2 (x))| 2 . 



(19) 



We now turn our attention to the third moments of A(x). 
There are two independent third moments (or triplets lfl6l[85l ) 
of A(x). In this article we only consider their diagonal ele- 
ments, and we define 

A(x) = <A*(x)A(x)A(x)>, (17) 

and 

X (x) = <A(x)A(x)A(x)>. (18) 

In Fig. [TJb) we plot the functions A(x) (solid line) and ^(x) 
(dashed line), evaluated on the z — plane and azimuthally 
averaged. Like the function k(x), both these moments are 
anomalous, but we find that they are real to within statistical 
uncertainty. Both functions reside primarily in the central re- 
gion where the condensate exists and, like p(x) and k(x), each 
exhibits its maximal absolute value at some distance from the 
origin. Although similar correlations have been discussed in 
the literature previously lfT6l l23 25 , 85, 86], we are not aware 
of any calculated forms for these functions with which to com- 
pare the results presented here. 

We stress that the expectation values A(x) and ^(x) vanish 
implicitly in the self-consistent mean-field (HFB and HFB- 
Popov) theories, and their appearance here is a signature of 



The differences between the actual fourth moments and their 
approximate Gaussian factorizations are quantified by the 
fourth cumulants of the field ll9"T| . There are three indepen- 
dent fourth cumulants of A(x), corresponding to the three in- 
dependent fourth moments of A(x). We consider here only the 
diagonal elements of these cumulants: 

<|A(x)| 4 > c = <|A(x)| 4 ) - 2<|A(x)| 2 ) 2 - |<A 2 (x)>| 2 , (20) 
<A*(x)A 3 (x)> c = <A*(x)A 3 (x)> - 3<|A(x)| 2 ><A 2 (x)>, (21) 
<A 4 (x)) c = <A 4 (x)> - 3<A 2 (x)) 2 . (22) 

In Fig. [TJc) we plot the azimuthally averaged values of 
<|A(x)| 4 > c (solid line), (A*(x)A 3 (x)> c (dashed line), and 
(A 4 (x)) c (dot-dashed line) on the z = plane. We note 
that these cumulants are small compared to the corresponding 
fourth moments. For example, from Fig.[T|a) it can be inferred 
that the maximal value of the moment (|A(x)| 4 ) is ~ 10~ 5 r~ 6 ; 
i.e., ~ 10 times the maximal value of (|A(x)| 4 ) c . We note 
that the fourth cumulants are somewhat more "noisy" than 
the lower correlations of the field; it is of course intuitively 
clear that the statistical demands for the accurate evaluation 
of cumulants will increase with the order of the cumulant. 

The cumulants presented here do, however, indicate that the 
fluctuations of A(x) exhibit a clear deviation from Gaussian- 
ity. Of particular interest is the cumulant (|A(x)| 4 ) c , as this 
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measures the deviation of the fourth moment (|A(x)| 4 ) from its 
Wick factorization [Eq. tfl9j]. This factorization becomes ex- 
act in the limit that the noncondensate fluctuations are Gaus- 
sian, and motivates the approximate factorizations of field- 
operator products made in mean-field theories in order to re- 
duce the second-quantized Hamiltonian to a (self-consistent) 
quadratic form [l4|j. In fact, Morgan ||20ll has shown that the 
treatment of the Hamiltonian term quartic in field operators 
in his second-order self-consistent calculation is equivalent to 
the Wick factorization employed in the HFB theories. Effects 
which lead to the appearance of a nonzero fourth cumulant 
(|A(x)| 4 ) c at equilibrium thus enter at a higher order of per- 
turbation theory. We note that the cumulant (|A(x)| 4 ) c found 
here is negative, which can be understood intuitively: the field 
admits non-Gaussianity of its fluctuations so as to lower the 
interaction energy associated with a product of four noncon- 
densate operators below the level obtained in self-consistent 
mean-field approaches. 



B. Many-body interaction effects 

We now relate the anomalous and non-Gaussian correla- 
tions of the classical field discussed in Sec. 



IV A 



to many- 
body interaction effects in the field which provide corrections 
to the mean-field potentials experienced by the condensate 
and noncondensed atoms. 



I. Mean-field potentials 

It is well known that the appearance of the anomalous av- 
erage in the generalized Gross-Pitaevskii equation of the HFB 
theory introduces many-body interaction effects which mod- 
ify the collisions between condensate atoms [20, 22, 23, 25 j. 
The origin of this effect is the Bogoliubov pair-excitation pro- 



cess (Sec. IV A i, which introduces the possibility that two col- 
liding condensate atoms are both scattered into (possibly oc- 
cupied) noncondensate modes, whereafter they may scatter off 
one another an arbitrary number of times before returning to 
the condensate. In general, the entire series of such "ladder 
diagram" processes 071 contributes to the effective interac- 
tion between condensate atoms in the finite-temperature sys- 
tem. In the HFB theory lTl4"l . it is found that the condensate 
<t>o obeys a generalized Gross-Pitaevskii equation of the form 
[Hs V -ii+Ua\<& {) \ 2 +2UQp\Q>Q+U()K<&l = 0, and indeed this form 
is also obtained in the formalism of Ref. EOll . The analyses 
of Refs. lF20ll25l show that the term Uok<$>q corresponds to the 
introduction of the many-body T matrix in the description of 
condensate-condensate interactions. We therefore reformulate 
this term as an additional potential experienced by the con- 
densate, which then obeys the nonlinear eigenvalue relation 
£ HFB d) = [i/ sp - / /+f/ l<l>ol 2 +2C/oP + f/o( l I>pVl<l>ol 2 ]<I) () = 0. 
The mean-field potential experienced by the condensate due to 
its own self-interaction 
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FIG. 2. (Color online) Many-body interaction effects in the classical- 
field equilibrium: (a-b) Mean-field potentials experienced by the 
condensate and their many-body corrections, (c) Local chemical po- 
tential of the classical-field condensate, in varying degrees of approx- 
imation (see text), (d) Effective interaction strengths deduced from 
the correlations of A(x). Gray solid lines in (b) and (c) indicate the 
spatial density profile of the condensate. 



thus receives a correction 



AV c (x) 



U 



|4>o(x) 



r Re{[(DS(x)] 2 <A(x)A(x))}. (24) 



In Fig. |2ja) we plot the bare mean-field potential V® 
(blue/black solid line), the many-body correction AV C (dashed 
line), and their sum, the corrected condensate mean-field po- 
tential V c = V c ° + AV C (dot-dashed line). We indeed find that 
the presence of the anomalous average corresponds to a no- 
ticeable reduction of the mean-field potential £/ol < E»ol 2 - 

Similarly to the mutual interaction of condensate atoms, 
the scattering of condensate atoms by noncondensate atoms 
is also affected by many-body processes in the finite- 
temperature Bose field: Upon its interaction with a non- 
condensed atom, a condensate atom may be scattered out 
of the condensate, and undergo an arbitrary number of (re- 
peated) interactions with the noncondensate atom before re- 
turning to the condensate. The consideration of all such 
processes amounts to the inclusion of all ladder diagrams 
in the condensate-noncondensate interaction. In an analy- 
sis beyond the usual mean-field (HFB) approach, Proukakis 
et al. have shown iTZBI that the condensate obeys an equa- 
tion of motion of the form id,<f>o - -Chfb^o + UoA(x, f), 
where A(x, t) is the (time-dependent) triplet correlator corre- 
sponding to Eq. (17 1. The adiabatic elimination of A(x, t) 



V?(X) = f/o|3>o(x)| 2 , 



(23) 



from this equation of motion leads to the introduction of the 
many-body r -matrix in the description of the condensate- 
noncondensate interactions, similar to the role of k(x) in the 
generalized GPE of the HFB theory. The term UqA can there- 
fore be reformulated as an additional mean-field potential ex- 
perienced by the condensate, yielding the equation of motion 
id,<3>o = [£hfb + f/o ( l ) o^/l < l ) ol 2 ] ( l ) o- The bare mean-field po- 
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tential experienced by condensate atoms, due to the presence 
of the noncondensate, 



V n ° c (x) = 2f/ <|A(x)| 2 >, 
thus receives a many-body correction given by 11251 



(25) 



AV nc (x) = 



U 



<!'()( x) 



-Re{<£;(x)<A*(x)A(x)A(x))}. (26) 



It is important to note that while the correction Eq. ( 24 1 arises 
due to the Bogoliubov processes which are accounted for to 
all orders by the Bogoliubov particle-hole structure of excita- 
tions in the HFB theories [20 1, the correction of Eq. ( |26] > cor- 
responds to processes in which atoms are exchanged between 
the condensate and its excitations, which are not described 
in such theories. In Fig. |2jb) we plot the bare potential V^ c 
(blue/black solid line), the correction AV nc (dashed line) and 
the corrected potential V nc = V^ c + AV nc (dot-dashed line). 
We find that the correction is localized near the boundary of 
the condensate, where the noncondensate density is maximal, 
however, it is nonzero at smaller radii also. Again, we see the 
effect of the many-body correction is to weaken the mean-field 
potential experienced by the condensate. 



2. Local chemical potential 

The properties of a superfluid are intimately connected to 
the phase of the order parameter (i.e., the condensate), as the 
superfluid velocity is proportional to the gradient of the order- 
parameter phase: v s = (h/m)V9 f75l [93l [94). The chemi- 
cal potential of the superfluid is, in general, a spatially vary- 
ing quantity, corresponding to (ft times) the phase-rotation 
frequency of the condensate |95l . A gradient in the chemi- 
cal potential therefore corresponds to an acceleration of the 
superflow, and a stationary superfluid should exhibit a spa- 
tially uniform chemical potential ||75l . Here we define a lo- 
cal condensate chemical potential corresponding to the ac- 
tion of the aforementioned (generalized) mean-field operator 
-Chfb + t^o^o^/l^ol 2 on the Penrose-Onsager condensate: 

/ic(x) =[// sp <Do(x)]/<D (x) + t/o[|3>o(x)| 2 + 2<|A(x)| 2 >] 

1 "° ? Re{[(D;(x)] 2 <A 2 (x)) + <DS(x)<A*(x)A 2 (x)>}. 

(27) 



l<5 (x)| 2 



We note that this local chemical potential is analogous to that 
derived by Zaremba et al. (87 1 within a symmetry-breaking 
framework. By evaluating p c (x) and approximations to it that 



neglect the many -body correction terms appearing in Eq. ( 27 1, 
we can assess the importance of these terms. Specifically, in 
addition to the full form of yU c (x), we consider the chemical 
potential obtained by evaluating only those quantities on the 



first line of Eq. ( 27 1 [which we will refer to as the "Popov" 
chemical potential, yUp op (x)], and that obtained by evaluating 
all but the last term of Eq. ( 27 1 [the "HFB" chemical potential, 
Mhfb(x)]. 



We calculate the action of H sp on the condensate <&q in the 
spectral representation, and azimuthally average the resulting 
quantity [H sp Oo(x)], and the individual factors of the remain- 
ing terms on the RHS of Eq. ( 27 1, on the z — plane in order 



to improve the estimates of these quantities and obtain radial 
representations of yu c (x), yU Pop (x), and jU H fb(x). The results are 
presented in Fig. |2jc), along with the effective eigenvalue Aq 
of the condensate orbital, obtained from the temporal analysis 
of Ref. 1 54 1 . We note first that the Popov-level chemical po- 
tential yKp op (x) (thick solid black line) is somewhat nonuniform 
and overestimates the condensate eigenvalue Aq (thin solid 
black line) over the extent of the condensate mode [where the 
concept of /U c (x) is meaningful]. The HFB chemical potential 
A*hfb(x) (dashed line) is smaller than yUp op (x) throughout the 
extent of the condensate mode, and approaches Ac much more 
closely at the smallest radii, while increasingly overestimating 
Aq at increasing radii [96 1. 

By contrast, the full chemical potential p c (x) (dot-dashed 
line) is much more uniform and in good agreement with 
the eigenvalue Aq aside from small discrepancies at the in- 
ner and outer extremes of the condensate orbital. This result 
shows that the condensate orbital obtained from the Penrose- 
Onsager analysis is consistent with a mean-field picture of 
the condensate and the complementary thermal component 
of the field, provided that the corrections to the mean-field 
potentials due to many-body effects are taken into account. 
This explains the essentially uniform phase rotation of the 
(Penrose-Onsager) condensate mode discussed in Ref. If54l . 
and moreover demonstrates the consistency of the equilibrium 
classical-field formalism with other formulations of the finite- 
temperature Bose-gas problem (e.g., Ref. Il87l ). We stress 
that the quantities jUp op (x) and /Uhfb(x) do not correspond to 
the condensate chemical potentials that would be obtained 
in formal calculations using the corresponding self-consistent 
mean-field formalisms; such calculations would of course de- 
scribe a condensate mode with a uniform chemical potential. 
Nevertheless, the results presented here give an indication of 
the order of the error involved in calculations employing such 
approaches. 



3. Effective interactions 

We now reconsider the corrections to the mean-field poten- 
tials experienced by the condensate (and noncondensate) in 
terms of effective interaction strengths, which feature in the 
theory introduced in Refs. Il22ll97ll . and reviewed in Ref. 11231 . 
Briefly, this so-called "gapless-HFB" (GHFB) theory is based 
on the identification that the inconsistency of the HFB theory 
arises because the theory introduces many-body corrections 
to the condensate-condensate scattering, but that effects of 
this order (i.e., the many-body T-matrix approximation |[T8l ) 
are not included in the description of interactions between the 
condensate and its excitations. The GHFB theory (of which 
there are two variants [22 1) thus introduces spatially depen- 
dent effective interaction strengths into the HFB -Popov theory, 
so as to approximate the effects of the many-body T-matrix 
corrections. 
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In Fig. |2jd) we plot three effective interaction potentials de- 
duced from the results of our classical-field simulation. We 
define a condensate-condensate effective interaction by group- 
ing the condensate-condensate mean-field potential V<? with 
its many -body correction AV C , and absorbing the effect of the 
latter into a spatially dependent redefinition of coupling con- 
stant [23: 



C/ C _ C (X) = Uq 

Making the substitution Uq — 
mean-field potential V® [Eq. 
the corrected potential V c 
action t/ c _ c (x) plotted in Fig 
an alternative visualization of 
potential presented in Fig.|2Ja 
with the effective interactions 



1+ (AA>\ 



d> 2 



(28) 



t/ c _ c (x) in the bare condensate 



(23 i] we obtain (by definition) 
Vf + AV C . The effective inter- 
|2Jd) (solid line) thus provides 
the correction to the mean-field 
), and is qualitatively consistent 
presented in Refs. 11221 |23l l84l . 



to noncondensate atoms, we identify 



U ac - ac (x) = U \l + 



<IM 4 )c 
2<|Ap>2 



(30) 



We observe that this quantity [dot-dashed line in Fig. |2Jd)], 
like U c - C and £/ c _ nc , describes interactions which are sup- 
pressed most strongly around the condensate periphery. How- 
ever, the maximal level of suppression here is of the order 
of ~ 5%, smaller again than that in the case of condensate- 
noncondensate scattering. The small value of this correction 



is to be expected, as Eq. ( 30 1 represents an average over all 
(coherent-region) excited-state collisions, many of which take 
place at high energies, where the many-body T matrix reduces 
to the two-body T matrix 112011221 . 



V. DEPENDENCE OF CORRELATIONS AND 
INTERACTION POTENTIALS ON THE FIELD ENERGY 



Next, we define an effective interaction appropriate to 
condensate-noncondensate scattering by absorbing the many- 
body correction AV nc to the mean-field potential V" c experi- 
enced by the condensate, due to the presence of the noncon- 
densate, into a local coupling constant 



In this section, we consider the dependence of the non- 
condensate correlation functions and interaction potentials de- 



t/ c - nc ( X ) EE Uq\ 1 + 



(A*AA) 
2O <A*A> 



(29) 



Making the substitution Uq -> U c - nc (x) in V° c [Eq. (25 i] 
clearly yields the corrected mean-field potential V nc = VjJT + 
AV nc - We observe that f/ c _ nc [dashed line in Fig. [2jd)], like 
I/ C _ c , is suppressed below the uniform value of the bare in- 
teraction potential Uq, and that this suppression is most pro- 
nounced around r w 4 (i.e., the location of the edge of the con- 
densate and the peak of the thermal-cloud density). However, 
the suppression exhibited by t/ c _ nc is significantly less than 
that of U c - C , suggesting that the so-called G2 variant of the 
GHFB theory ||22||23]|27), in which U is replaced by C/ C _ c ev- 
erywhere it appears in the Gross-Pitaevskii and Bogoliubov- 
de Gennes equations, significantly overestimates the suppres- 
sion of the noncondensate mean-field potential experienced 
by the condensate. By contrast, the Gl variant of the theory 
retains the unmodified coupling constant Uq in the mean-field 
potential V nc , and thus neglects any suppression of this poten- 
tial. 



Finally in this section, we introduce an effective interac- 
tion strength for noncondensate-noncondensate interactions. 
To define this quantity, we recall that the appearance of a neg- 
ative fourth cumulant (|A| 4 ) C implies that the noncondensate- 
noncondensate interaction energy £4 = (Uq/2) Jdx{\A\ 4 ) is 
suppressed somewhat below the value obtained by the Gaus- 
sian factorization of the moment (|A| 4 ) assumed in the HFB 



theory (Sec. IV A 1. By assuming that the difference between 
£4 and the corresponding approximate Gaussian expression 
can be accounted for by an effective (spatially-dependent) in- 
teraction strength associated with the mean-field potential due 



scribed in Sec. IV on the total energy (and thus temperature) 
of the classical-field equilibrium. 



A. Correlation functions 

In Fig. [3] we plot the spatially integrated values of 
correlation functions of the noncondensate field, obtained 
from PGPE equilibria with energies in the range E e 
[9.5,24.0]N c Hu> r . The temperature and chemical potential of 
the field, obtained using the microcanonical thermodynamic 
formalism of Rugh i55l l67l l70l . are presented in Fig. [5|a) 
for reference. In the inset to Fig. |3ja), we plot the quantity 
Cb = (|Q'oI 4 )/(|Q'oI 2 ) 2 as a function of the field energy. Bezett 
and Blakie [80 1 have suggested this quantity as the appropri- 
ate generalization of the Binder cumulant associated with the 
phase transition of the homogeneous gas ll67ll to the harmon- 
ically trapped case. Campostrini et cd. [98| found the critical 
value (CbW = 1.243 at the phase transition of the classi- 
cal three-dimensional XY model, and this value was used to 
identify the critical temperature in PGPE simulations of the 
homogeneous Bose gas in Ref. [67|. The solid line in the in- 
set to Fig. |3ja) interpolates between the values of Cb obtained 
from our simulations (circles), and intersects the critical value 
(CB)crit (horizontal dashed line) at E — 21.1N c ha> r (vertical 
dashed line), which we take as an estimate of the critical field 
energy for this system. The critical energy is indicated by a 
vertical dashed line in Fig.[3|a-d). 

Figure |3jb) shows the spatially integrated normal ther- 
mal density Jdxp(x) (circles) and anomalous thermal den- 
sity fdxi<(x) (squares). The results are similar to those 
obtained in ll54l . in which the anomalous correlation func- 
tion k(x) was defined in terms of a symmetry-breaking in- 
terpretation of the classical-field trajectories. It is impor- 
tant to note that the noncondensate population increases ap- 
proximately linearly with energy (and thus temperature) up 
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where critical fluctuations subside and modes recover Gaus- 
sian fluctuations, A(x) vanishes. 

In Fig. |3jd) we plot the integrated values of the fourth 
cumulants of the noncondensate field A(x). The cumulant 
(A 4 (x)) c behaves similarly to k(x), as expected, given that this 
cumulant quantifies the difference between the purely anoma- 
lous (condensate-induced) moment (A 4 (x)) and its naive fac- 
torization in terms of k(x) [see Eq. (|22]i]. Similarly to /t(x), the 
integrated value of (A*AAA) C vanishes at a slightly higher 
temperature. From the definition [Eq. (21 1] of this cumu- 
lant, it is clear that this behavior reflects that of the moment 
(A* AAA), which can be understood in a similar manner to 
the behavior of A(x). By contrast, the "normal" cumulant 
(|A| 4 ) C reaches its maximum absolute value around the tran- 
sition to the normal phase. This shows that many-body inter- 
action effects are most important in the noncondensate field 
in this regime, in which the condensate is surrounded by a 
much larger region of the field which exhibits suppressed den- 
sity fluctuations but no phase coherence (i.e., a quasiconden- 
sate 111 0011 ). as we discuss further in Sec.|VI| 



B. Interaction energies 



FIG. 3. (Color online) Dependence of thermodynamic parameters 
and noncondensate correlations on the classical-field energy, (a) 
Temperature and chemical potential of the field. Magnitudes of spa- 
tially integrated (b) second moments, (c) third moments, and (d) 
fourth cumulants of A(x). Inset: Dependence of the Binder cumulant 
Cb on the field energy. The solid line interpolates smoothly between 
the numerical data points (circles). The dashed lines indicate the crit- 
ical value (CB)crit = 1.243 and the corresponding critical field energy 
(see text). 



to the transition of the field to the noncondensed phase in 
the idealized "PGPE system" consisting of a fixed field pop- 
ulation distributed over the modally finite coherent region 
ll54l l55l l64l [671 . Upon the consideration of the above- 
cutoff fraction of atoms, one obtains the expected, geometry- 
dependent scaling of the (non-)condensate fraction with tem- 
perature (see, for example, Ref. Il53l ). By contrast, the inte- 
grated anomalous density reaches its maximum absolute value 
at intermediate energies (temperatures), and vanishes as the 
system approaches the phase transition, in agreement with 
previous studies 112211231 I54ll84t 

The integrated third moments J dx A(x) and J dxx(x) are 
plotted in Fig. [^c). The quantity Jdxx(x) reaches its 
maximum at intermediate energies, and vanishes along with 
J dxi<(x) at E » 2l.5N c ha> r . The quantity JdxA(x), however, 
remains nonzero at the transition, and tends to zero at a some- 
what higher energy. We reiterate that although the definition 
[Eq. ( [TT} ] of A(x) is singular when ao = 0, this condition 
does not occur in practical simulations |99|. Nonzero values 
for J dx A(x) are obtained at the phase transition as A(x) repre- 
sents the scattering of atoms into and out of the "condensate" 
mode ^o( x ). which persists in the critical regime, even though 
the mode;fo(x) is no longer coherent. At higher temperatures, 



We now consider the variation of the many -body interaction 



effects identified in Sec. IV B 3 as functions of the field energy. 
We introduce three interaction energies: 



E c . 



If 

!■ 



fl fxf/ c _ c (x)|O (x)| 4 , (31) 



<ixt/ c _ nc (x)|(Do(x)| 2 p(x), (32) 



t/ nc _ nc (x)p 2 (x), (33) 



corresponding to the condensate-condensate, condensate- 
noncondensate and noncondensate-noncondensate interac- 
tions, respectively. We compare these three energies with the 



corresponding uncorrected energies, obtained from Eqs. (31 



( 33 1 by replacing the spatially dependent interaction strengths 
[J/ c -c(x), etc] with the bare interaction strength Uo, which we 
denote by E®_ c , etc. 

In Fig. HI we plot the corrected interaction energies 
Eqs. (f3T]> - (f33j> as fractions of the corresponding uncorrected 
energies. The approximate location of the phase transition 
determined from the analysis of the Binder cumulant Cb 
(Sec. VAi is again indicated by a vertical dashed line. The 



absolute magnitudes of the corrected energies are shown in 
the inset to Fig. [4] Considering the condensate-condensate 
interaction energy E c - C (circles in Fig. |4j, we note that it con- 
tinues to decrease (relative to E®_ c ) as E increases towards the 
transition, and reaches its minimum at E » 2lN c haj r . We 
contrast this with the behavior of J dx k(x) [Fig. 3(b)], which 
reaches its maximum (absolute) value at £ a; 15NJtto,-. At en- 
ergies E > 23N c fi(L) r , where no condensate is present [i.e., the 



mode^o(x) is completely incoherent], the expression Eq. (31 
is meaningless, and so we do not present data for £ c _ c at these 
highest energies. We find that the total condensate-condensate 
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FIG. 4. (Color online) Dependence of the corrected interaction en- 
ergies on the energy of the classical-field equilibrium. Interaction 
energies are plotted as fractions of uncorrected interaction energies 
(see text). The absolute values of the interaction energies are shown 
in the inset. 



interaction energy is suppressed by ~ 25% at the transition, 
which is consistent with previous mean-field studies of har- 
monically trapped condensates 11221 l23l [84*1 . and should be 
contrasted with the complete vanishing of the condensate- 
condensate interaction at the transition temperature in the ho- 
mogeneous case l4l [T8l l26l . 

Turning our attention to the energy E c - nc of the condensate- 
noncondensate interaction (triangles in Fig. 0), we observe 
that Zic-nc /E°_ nc , like E c ^ c /E^_ c , reaches its minimum (~ 
0.88) around the transition to the normal phase. At higher en- 
ergies (temperatures), the interaction returns to its bare value; 
in this limit, the "condensate" mode xo( x ) is simply a ther- 
mal (Gaussian) mode (i.e., Cb 



2), and Eq. (32i reduces 
to the Hartree-Fock interaction energy between this mode 
and the rest of the field. The energy of the noncondensate- 
noncondensate interaction (squares) is similarly suppressed 
most strongly at around E * 2lN c ha> r , but is subject to less 
suppression than the other two interaction energies. 

We note the reasonable agreement between the critical en- 
ergy estimated from the consideration of the Binder cumulant, 
and the energy at which maximal suppression of the field in- 
teractions occurs, which provides an independent estimate of 
the critical point (see Ref. [22] and references therein). This 
suggests that the critical value (CB)crit = 1-243 appropriate to 
the phase transition of the homogeneous Bose gas does indeed 
yield a good estimate of the transition energy (temperature) in 
the harmonically trapped system. This is in contrast to the 
results of Ref. Il80l which suggested, based on an analysis of 
the field correlation length, that the critical value of Cb in the 
harmonically trapped system may be significantly lower than 
that in the homogeneous case. [Indeed the results presented in 
Figs.[3]and[4]of this work suggest that (CB)crit slightly under- 
estimates the correct critical value]. 

Finally in this section, we note that the interaction ener- 
gies considered here constitute an increasingly small fraction 



of the total field energy as the field energy (temperature) in- 
creases (inset to Fig. [4]). In particular, in the critical region, the 
noncondensate-noncondensate interaction energy is < 10% of 
the total field energy, and the remainder is mostly comprised 
of the kinetic energy of the noncondensate fraction. At the 
transition, the contribution of beyond-Gaussian fluctuations 
to the total field energy is therefore only of order ~ 1% (cf. 
the agreement of HFB-Popov and PGPE calculations of the 
critical temperature to about this level reported in Ref. [53 1). 



VI. SUPPRESSION OF DENSITY FLUCTUATIONS IN THE 

FIELD 



We have seen in Sections |IV] and [V] that the classical-field 
equilibria exhibit a negative fourth cumulant (|A| 4 ) C , corre- 
sponding to the suppression of density fluctuations in the non- 
condensate field A(x) below the Gaussian level assumed in 
the HFB and HFB-Popov formalisms. We now consider the 
suppression of density fluctuations in the total classical field 
(/'(x), and identify and compare the various contributions to 
this suppression. We follow Ref. [ 100 1 in considering the cor- 
relator 



Q(x) = 2<|,A(x)| 2 ) 2 - <|<A(x)| 4 >, 



and defining the quasicondensate density 



(34) 



« Q (x) = VGW = [2 - g (2) (x)] 3 <|«A(x)| 2 ), (35) 

where the classical coherence function g (n) = (\iff\ 2 ")/(\iff\ 2 )" 
fi4l 11011 [102 1. The quasicondensate density «q(x) has 
been used by Bisset et al. ll60l to characterize the super- 
fluid transition of the quasi-two-dimensional Bose gas (see 
also Ref. 1 103 1). We recall that the local correlation func- 
tion adopts values of = n! and = 1 for "nor- 
mal" thermal fields (those without anomalous correlations) 
and purely coherent fields, respectively 1144111041 . In Fig. [5ja) 
we plot the (azimuthally averaged) quasicondensate density 
«q for the PGPE equilibrium with energy E = 18.0N c n(L> r 
(blue/black solid line), along with the condensate density |<t>o| 2 
(gray solid line) and total field density (|iA| 2 ) (dashed line). 
We observe that the quasicondensate density Hq is everywhere 
greater than the condensate density, and in particular, «q re- 
mains significant at radii r > 4r , for which the conden- 
sate density essentially vanishes. In Fig.|5jb) we plot the in- 
tegrated condensate density, i.e., the condensate population 
«o = Jdx \<S>o\ 2 (circles), and the integrated quasicondensate 
density Nq = JdxriQ (squares), as functions of the classical- 
field energy E. We find that the quasicondensate population is 
greater than the condensate population for all energies consid- 
ered, and that the difference between the two populations in- 
creases with increasing E. Most notably, in the highest energy 
equilibria considered (E > 22.QN c Ha> r ), in which no conden- 
sation remains, the quasicondensate population is still signif- 
icant. The presence of this quasicondensate at high energies 
(i.e., in the absence of condensation) clearly corresponds to 



12 




1.5 3 4.5 6 
Radius (units of ro) 



10 12 14 16 18 20 22 24 
Energy (units of N^hus,-) 





Q 




(2-Cb)|4>„| 4 




2(*!) 2 (A a ) 




i*:(A*A J > 















1.5 3 4.5 6 
Radius (units of ro) 



Q 

„ (2-C B )|4 | 4 
4 -2(*;) ! {A ! ) 
* -4*;{A*A 2 ) 
■ -(|A|'>. 



10 12 14 16 18 20 22 24 
Energy (units of N^huj,-) 



FIG. 5. (Color online) Suppression of density fluctuations in the 
classical field, (a) Condensate, quasicondensate, and total density 
profiles in a representative classical-field equilibrium with energy 
E = l8.0N c Tia) r . (b) Dependence of the total condensate and qua- 
sicondensate populations on the field energy, (c) Dominant con- 
tributions to the quasicondensation correlator Q in the case E = 
\%.QNJioj r . (d) Dependence of the dominant contributions to Q on 
the field energy. 



the significant fourth cumulant (|A| 4 ) C found in this regime 
(Sec.|VAli. 

By expressing (\ift\ 2 ) and (\ift\ 4 ) in terms of aoXo and 5, sub- 
stituting 6 — > A and assuming that the condensate can be fac- 
tored out of these moments [e.g. <|o , oA'ol 2 |5| 2 ) ~ I^oRAI 2 ) 
and ((a^x* Q ) 2 6 2 ) ~ (<E>*,) 2 (A 2 )], we obtain the approximate 
form |105] 



|<AA>| 2 + 2Re{(€g 2 <AA>} 



-4Re{OS<A*AA» + <|A| 4 >c , (36) 



Q * (2 - C B )|O 



where the Binder cumulant Cb accounts for the amplitude 
fluctuations of the condensate mode (see Sec. |V A| >. We note 
that only by neglecting all but the first term on the RHS of 
Eq. (B6]l do we obtain (in the limit Cb — > 1) «q = | *3>o I 2 (cf. 
Ref. 1106]). The second term on the RHS of Eq. (36 1 repre- 
sents the enhancement of density fluctuations in the noncon- 
densate component relative to the normal Gaussian [g (2) = 2] 
level, due to the anomalous nature of its fluctuations (consider 
that <|A| 4 > * 2<|A| 2 ) 2 + |<A 2 >| 2 > 2<|A| 2 ) 2 ) HUT] fT08l . and 
thus subtracts from Q. In practice we find that the contribu- 
tion of this term is small. The remaining terms on the RHS of 



Eq. ( 36 1 all represent reductions in the density fluctuations of 
the field ^(x), ar, d thus a dd to the quasicondensate density. 

In Fig. [5jc) we plot the (azimuthally averaged) quantity Q 
for the case E - 18.0N c nu> r (solid blue/black line), along with 
the three dominant contributions to this quantity: (2-Cb)| ( I ) oI 4 
(solid gray line), -2(<D*) 2 <A 2 ) (dashed line) and -4<&;<A*AA) 



(dash-dot line). The contributions of |(AA)| 2 and (|A| 4 ) C are 
essentially negligible for this field energy, and the suppression 
of density fluctuations here is therefore primarily associated 
with the presence of a condensate in the field. However, the 
extent to which density fluctuations are suppressed is some- 
what greater than that due to the absence of density fluctua- 
tions in the condensate itself. This suggests that an estimation 
of the condensate density based on an experimental measure- 
ment of g (2) could overestimate the size of the condensate con- 
siderably. We plot the spatially integrated values of the corre- 
lator Q and the dominant contributions to it as functions of the 
field energy in Fig. |5Jd). We observe that the contribution of 
-2(<t>p) 2 (AA) (triangles) is the dominant beyond-condensate 
contribution to Q at the lowest energies, but reaches its max- 
imum magnitude at E as l3.5N c ftaj n and is overtaken by the 
contribution of -4<]>* ) (A* AA) (pluses) at higher energies. This 
contribution is itself overtaken soon after by the contribution 
of the fourth cumulant (|A| 4 ) C (crosses), which is of course the 
dominant contribution in the critical regime. 

The results presented here should be compared with the 
behavior of the homogeneous Bose gas in two dimensions, 
for which it is found that the quasicondensate density is in 
general greater than the superfluid density, and that the qua- 
sicondensate correlations appear at temperatures well above 
the critical temperature for the appearance of superfluidity in 
the system [100 10911 . Classical-field calculations for the har- 
monically trapped quasi-two-dimensional Bose gas appear to 
be consistent with a local-density version of this picture ll60l . 
It is perhaps not surprising that we observe similar behavior 
in the three-dimensional case considered here: the compara- 
tively large transverse susceptibility of the field [1 10] leads to 
the (local) destruction of order by phase fluctuations before 
size fluctuations become significant 1 1 11 1. 



VII. CONCLUSIONS 

We have presented a methodology for the calculation of 
general spatial correlations of the finite-temperature Bose 
gas, including anomalous correlations, in a classical-field ap- 
proach. We introduced U(l)-symmetric classical-field vari- 
ables, analogous to number-conserving quantum ladder op- 
erators considered by other authors P71 - I491 . in order to rig- 
orously quantify the anomalous correlations in the micro- 
canonical PGPE ensemble. We have demonstrated that the 
finite-temperature Bose field exhibits non-Gaussian correla- 
tions which are classical in origin, and by relating these cor- 
relations to many-body corrections to the mean-field interac- 
tion potentials C20l l22l 1251 1571 , we have explicitly demon- 
strated that the classical-field theory provides an intrinsic 
nonperturbative description of many-body processes in the 
field. We discussed the role of these processes in determin- 
ing the condensate mode, and demonstrated the consistency 
of the Penrose-Onsager 11721 definition of condensation in the 
classical-field equilibrium. Finally, we elucidated the contri- 
bution of the anomalous and non-Gaussian correlations of the 
noncondensate to the overall suppression of density fluctua- 
tions in the field, and discussed the distinction between the 
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quasicondensate defined by the suppression of density fluctu- 
ations, and the true condensate defined by the suppression of 
phase fluctuations. 

Our demonstration of the intrinsic description of many- 
body interaction effects provided by the classical-field method 
underlines the utility of these techniques in describing low- 
dimensional systems [60, 1 12], in which such effects can have 
profound consequences for the structure of the field equilib- 
rium II113H1151 , Indeed our results illustrate the complexities 
introduced by the independence of phase and density fluc- 
tuations Bllll II 161 in the inhomogeneous finite-temperature 
Bose gas even in the comparatively straightforward three- 
dimensional case. Moreover, although we have only consid- 
ered systems at equilibrium, our results demonstrate that the 
classical-field model of the low-energy coherent region does 
describe the higher-order processes associated with quasipar- 
ticle damping, and the exchange of atoms between the con- 



densate and its excitations, which are essential for the descrip- 
tion of nonequilibrium dynamics of the finite-temperature 
Bose field lfT2ll2Tll25l 1291. 
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